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Comment on "Dynamical Heterogeneities in a 
Supercooled Lennard-Jones Liquid" 

In two recent interesting letters Q evidence was pre- 
sented for the existence of a growing dynamic correla- 
tion length when we approach the glass transition from 
the liquid phase (a similar divergence is present also 
in the off-equilibrium dynamics below T c Here 
we would like to point out that this phenomenon can 
be easily predicted using the replica approach of ref- 
erence [pi. Let us call p(x,t) the density at time t 
(p(x,t) = Yli^( x i(t) ~ x ))- K is convenient to con- 
sider the quantity Q Zl (x,t) = p(x — z\/2, 0)p(x + Zi/2, t). 
Using it we can construct the correlation C{x — y,t) — 
{Qz-i (x, t)Q Z2 (y, £)) that should be weakly depend on the 
values of z\ and Z2 for \z\\, \z 2 \ << x — y. The data in 
jjj can be read as evidence for the existence of a correla- 
tion length that controls the decay of C at large spatial 
argument. This length reaches a maximum at a time t 
comparable with the relaxation time of the system and 
this time grows for lower temperatures. Let us take the 
mode coupling point of view, assuming an incipient diver- 
gence of the relaxation time and breaking of ergodicity 
at a temperature T c . Exactly at the mode coupling tran- 
sition we would sample at large times a single ergodic 
component and obtain 


(p{x- 


Coo (x) = limt^oo C(x —y,t) = 
Zl /2)p((y + z 2 /2)) c {p{x - Zl /2)p{y 


V2)c (1) 


where the expectation value is taken inside the same er- 
godic component and c stands for connected. As soon 
as we approach T c and the times become large, the dy- 
namic correlation function goes, in the appropriate time 
window, to the critical one. Our burden is to prove that 
there are (within the approximation used) long range cor- 
relations at the mode coupling point. 



FIG. 1. We plot dW/de as a function of the coupling e for 
the model in [§] and T = T c (l-0.l) (dashed line), T = T c (full 
line) and T = T c (l + 0.5) (dotted line). Notice the divergence 
of the d 2 W/de 2 for e -> for T = T c in the high branch of 
the curve. 


may be characterized as follows. Consider a system that 
thermalizes at temperature T with an attractive cou- 
pling e with a reference configuration which is itself a 
fixed equilibrium configuration at temperature T. We 
call W the free-energy as a function of e. There are three 
possibility for the expectation value of (dW/de) \ e =o- ( a ) 
A non zero stable value below the ideal glass transition 
(T <T K ) (b) Flatly a zero value above the Mode Cou- 
pling transition (T > T c ) (c) A metastable non zero value 
in the region Tk < T < T c . The metastable state ex- 
ists also for non zero coupling and is destabilized on the 
spinodal line. At T c the spinodal point is at e = 0. The 
situation is sketched in fig 1. It is a general feature of 
spinodal points that the second derivative of the free en- 
ergy, and its associated correlation length are divergent. 
This is precisely the length associated to Coo(x). The 
learned reader may notice that we are just finding what 
one should expect from the "principle of marginal stabil- 
ity" [[| which states that the mode coupling transition 
point is critical. Similar conclusions follows also from 
the geometrical description of [Q. It should also be ev- 
ident that in principle this correlation length could be 
computed directly in the mode coupling approach. Un- 
fortunately such a computation is not so simple and the 
method of || provides an useful short cut. We stress 
that the divergence is present only in the average of the 
product of two correlation functions inside the same er- 
godic component. No divergent length would be seen in 
equilibrium correlation functions, whose dominant con- 
tribution comes from configurations in different ergodic 
components. 
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In H it was shown that the mode coupling transition 
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